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<N . Abstract 

^\ n 

^H, A recent experiment [l[ suggested tliat bisniutii forms liexagonal close packed (HCP) films on 

the Ni(lll) surface, of heights 3, 5 and 7 layers. A quantum size effect based on free electrons was 



proposed in explanation. To test this idea, we calculate the total energies of Bi on the Ni(lll) 
surface using density functional theory. We find that HCP film stabilities disagree with the observed 
odd layer preferences, and the structures are mechanically destabilized by adding capping atoms 



C/3 

1^ \ which pucker the HCP layers. Furthermore, we find that rhombohedral films based on the bulk 



Bi structure are energetically more favorable than the proposed HCP films. These structures also 
favor odd numbers of layers, but owing to covalent chemical bonding rather than confinement 



'^ \ of free electrons. Specifically, a strongly bound adsorbed surface monolayer forms, followed by 

O ' bulk-like rhombohedral bilayers. 



PACS numbers: 



INTRODUCTION 

During the growth of thin metaUic films, confinement of electrons can favor film heights 
that are commensurate with half the Fermi wavelength [2]. This well known quantum 
size effect due to confinement (QSE) has long been studied in various materials, including 



bismuth. Bismuth differs from conventional OS 
especially long Fermi wavelength (~ 40 nm) |3|, 



elements because it is a semimetal with an 
J] that causes various physical properties of 



Bi films to oscillate with long periods [Sj. In ultrathin films (~ 1 nm), Bi exhibits allotropic 
transformations from puckered pseudocubic films to hexagonal symmetry films on Si(lll) 
and on some quasicrystal surfaces [6|, l7|. Experiments and first principles calculations concur 
that both film types exhibit bilayer growth, due to the exotic bonding character of Bi, rather 



B 



than QSE l8|. In contrast, the initial growth of Bi films on metallic substrates has not been 

n 

well studied until a recent experiment [JJ] reported stable 3, 5 and 7 layer Bi hexagonal close 
packed (HCP) films of Bi on the Ni(lll) surface. They attempted to explain the stability 
by quantum confinement based on a free electron model. A more in-depth theoretical study 
is needed to understand this possible short period QSE in Bi. 

Bulk Bi takes the rhombohedral structure of Pearson type hR2 (prototype a- As) common 
to group-V semimetals, which is distorted from the simple cubic structure by a Jones-Peierls 



Is 



mechanism (9|. The bulk Bi structure is best described as a stacking of bilayers in the [001] 
direction. Here we use 3-element hexagonal indexing, it would be [111] using rhombohedral 
indexing. Each bilayer has height 1.43 A and is separated from the adjacent bilayer by 
2.51 A. Within the (001) plane the Bi spacing is 4.53 A. However, three strong covalent 
bonds of length 2.98 A link each Bi atom to others within each bilayer, while three weak 
metallic bonds of length 3.62 A connect each Bi atom to others in the adjacent bilayer. The 
(001) plane is thus a natural cleavage plane, with divisions expected between bilayers. In 
addition to hR2, bulk Bi possesses many allotropes, especially at high pressure, including 
Pearson structure types mP4, mC4, cPl and cI2, but it does not take the HCP structure. 

n 

In their experiment pj], the authors grew Bi on a Ni(lll) surface at a temperature of 
473 K via vapor deposition. Based on the low energy electron diffraction (LEED) patterns, 
deposition and film growth rates, and low energy electron refiectivity spectra, the authors 
proposed that initially a (7x7) wetting layer forms (hereafter referred to as adsorbed surface 
monolayer), which transitions to a 3 layer HCP film with a (3x3) surface cell sitting directly 



on the Ni substrate, as well as a 5 layer HCP film with a [3-112] surface cell, as coverage 



grows. However, at 422 K, Bi formed a 7 layer film with an (8x8) cell 10| surrounded by 
the 3 layer (3x3) film. The proposed 3 layer and 5 layer Bi films on Ni have in-plane lattice 
constant 3.7-3.8 A (See Appendix Al. for discussion of lateral strain). Taking their measured 
lattice constants and an assumed free electron valence of 5, they calculated that their 3, 5 
and 7-layer films were, respectively, 2.5, 4.0 and 5.0 Fermi wavelengths in height. As these 
structures and bond lengths have not been previously observed, and Bi is notoriously not 
free electron-like, a first principles electronic structure investigation is warranted. 

METHODS 



We apply electronic density functional theory, using the Vienna ab-initio simulation pack- 
age (VASP I 111 . Il2l |) to solve the Kohn-Sham equations with the Perdew-Burke-Ernzerhof 



3,Q 



(PBE [13|) parameterization of the generalized gradient approximation (GGA) for the el 
tron exchange correlation potential. We use projector augmented wave potentials 
with a fixed energy cutoff of 269.5 eV (the default for Ni). The d semi-core levels of Bismuth 
are treated as valence electrons. CoUinear spin polarization is used since Ni is ferromag- 
netic, though we test noncollinear magnetism to check the importance of spin-orbit coupling 
(SOC) for some structures. In the noncollinear calculations, we take the relaxed structure 
from collinear calculations and perform a static calculation. 

We construct models based on four Ni layers normal to the (HI) surface with Bi films on 
one side. Our cells include 22 A of vacuum, with periodic boundary conditions. Electrostatic 
energy created by the asymmetric charge distribution in the presence of Bi is small relative 
to the differences of surface energies. All structures are relaxed holding the cell sizes and 
bottom layer Ni atoms fixed, with in-plane lattice constants set by the relaxed bulk Ni 
structure. Energy convergence is carefully checked with respect to the vacuum size, k-point 
mesh and the number of Ni layers. 

RESULTS AND DISCUSSION 

We first compare the relative surface energies of the proposed Bi HCP films of different 
thickness on Ni(lll). Several quantities are needed to define relative surface energy: the 
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FIG. 1: (Color online) Relative surface energies of HCP Bi films with (3x3) and [3-112] Ni (111) 
surface cells. Note the energy conversion factor is 1 eV/A^=16.0 J/m^. Red points are from 
collinear calculations, blue points are noncollinear. Spin-orbit coupling affects the total energy 
quantitatively but does not alter the relative film stabilities. Data points connected by line segments 
correspond to coverages of integer numbers of monolayers (1-9 for 3x3, 1-11 for [3-112]). Extra 
data points in the (3x3) cell are three layer films with one extra capping atom at a valley site. 
Extra points in the [3-112] cell are 4 layer plus one atom, and 8 layer plus one atom at valley sites. 

total energy Efot of the Nsi atoms of Bi on the surface of the 4-layer Ni slab; the slope 
Esi, which is the linear part of Efot as it depends on Nsf, the energy E^'^'' of the Ni slab 
including its two free surfaces, each of area A. The slope Ebi can be considered as the energy 
of bulk Bi in the HCP structure with the in-plane lattice constant determined by the surface 
cell 16|. With these definitions, relative surface energy is 

[Etot - E'^f - EB,NB^]/A. 



Jrel 



(1) 



Fig. [U shows the relative surface energies of (3x3) and [3-112] HCP films. A structure is 
relatively stable if the second derivative of its relative surface energy with respect to film 
thickness is positive [17|. Relative stability occurs for (3x3) cells of 1, 3, 4, 6 and 8 layers 
and for [3-112] cells of 1, 4, 6, 8 and 10 layers. These predicted stabilities disagree with the 
experimentally discovered odd layer preference. 

Also shown in Fig. [T]are data points for films with single Bi adatoms. These data points 
lie below the relative surface energies of integer layer coverages, revealing that terminating 
on complete HCP layers is unfavorable. Strong puckering of the HCP layers occurs during 
relaxation of structures with adatoms. Relative surface energies calculated with SOC are 



shown for the (3x3) films. SOC influences the relative surface energies quantitatively but 
does not alter the sequence of stable structures. 

The relatively stable 4-layer (3x3) structure is illustrated in Fig. |21 The first layer Bi 
atoms strongly bond with the Ni surface atoms. On top of the first layer, Bi atoms form 
slightly puckered layers in order to achieve short Bi bonds with adjacent layers. These short 
covalent bonds lower the relative surface energy. 

The asymptotic period 2 oscillations of relative surface energies are reminiscent of the 
QSE prediction, but the minima occur for even numbers of HOP layers rather than odd. 
Moreover, we have done a genuine first principles study of QSE in free standing HCP Bi 
films (see Appendix A2.) and find that the actual predicted oscillation period is close to 3 
layers, not 2. Thus our first principles calculations cast doubt on both the HCP structure 
model and the proposed explanation in terms of QSE. 
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FIG. 2: (Color online) Side view of relaxed relatively stable 4-layer HCP Bi film on Ni(lll) (3x3) 
cell (units are A). Chemical bonding shown at right (Slightly tilted). 



We now seek alternate film structures. Inspecting Fig. |2] we note the interaction of the 
upper layer Bi with the surface monolayer is weak, so strong deviations of structure and 
bonding from bulk Bi due to the Ni substrate are not anticipated. Thus bulk-like hR2 
Bi (001) films are good candidates to grow on top of the surface monolayer. Henceforth, 
when we say "hR2 film" we include the surface monolayer. For example, a 3 layer hR2 
film (see Fig. [3]) consists of the surface monolayer plus an hR2 (001) bilayer. The in-plane 



lattice constant of bulk hR2 Bi is 4.53 A(see Appendix Al.). Meanwhile, \^ times the 
Ni interatomic spacing is 4.31 A which differs by only 5%. Furthermore, in our calculated 
energy of strained free standing Bi (001) films, the energy minimum occurs at in plane lattice 
constants 4.3 A, 4.4 A and 4.5 A for one, two and three bilayer films respectively. Thus the 
stable thin bilayer films match with Ni very well. 
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FIG. 3: (Color online) Side view of relaxed Bi surface monolayer and rhombohedral bilayer on 
Ni(lll) (3x3) cell (units are A). Chemical bonding shown at right (Slightly tilted). 

To illustrate the relative stability of various structures, we compare the surface enthalpy 
of formation, which is defined as. 



AH /A = [E, 



tot 



Tpslab 



rpbulk 



NB^]/A 



(2) 



which differs from the relative surface energy '^^ei in E<1- E only in our choice of reference 
energy for pure Bi, -E^"^ is the relaxed Bi bulk energy in the hR2 structure. 

Fig. m shows the enthalpy of formation for various surface structures with different thick- 
ness. Notice that both the (8x8) and (7x7) surface monolayer structure touch the convex 
hall which implies they are both energetically stable. The stable (3x3) surface monolayer 
is equavilent to both 1 layer HCP film and hR2 film with (3x3) cell. Moreover, for higher 
coverage, hR2 films have much lower energy than the HCP films. For instance, at the same 
total coverage of 16/9=1.8 Bi/Ni, the total energy the five layer hR2 film is 1.5 eV lower than 
the four layer HCP film on the bare (3x3) surface. This energy difference is much larger 
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FIG. 4: (Color online) Enthalpy of formation. Dot-dashed green points are the (sub)monolayer 
structures with (3x3) cell. Orange star and plus are the (7x7) and (8x8) surface monolayer 
respectively. Red points are HCP films and black are hR2 Bi films on the (3x3) cell. Blue points 
are HCP Bi films on [3-112] ceh. 

than the thermal energy, A;sT=40 meV at 473 K. The hR2 films are thus much more likely 
to form than the HCP films. The hR2 films favor odd number of layers (surface monolayer 
-|- integer bilayers) which is consistent with the experimental observations of 3, 5 and 7 layer 
films. However, the stability is due to the exotic chemical bonding of Bi rather than QSE. 
To further illustrate the stable sequence from equilibrium thermodynamics, we calculate 
the surface free energy. This quantity is the Legendre transform of the enthalpy of forma- 
tion (Eq. [2]), replacing the surface coverage with relative chemical potential AfiBi- From 
equilibrium thermodynamics, the most stable structure at a certain Bi chemical potential 
AfiBi minimizes the surface free energy 18|, 



J=[AH-A^lB^NBi]/A 



(3) 



where AH is the enthalpy of formation (Eq. [2]), AfiBi = f^Bi — E^m^ is the Bi relative 
chemical potential. The results in Fig. [5] shows the stable sequence is from bare surface to 
one atom on the (3x3) cell at AyUBi=-l-16 eV, to three atoms on the (3x3) cell at -1.11 eV, 
to (8x8) surface monolayer at -0.83 eV to four atoms ((3x3) surface monolayer) at -0.82 eV, 
to the (7x7) surface monolayer at .01 eV, and then finally to the infinite height hR2 films 
at around 0.06 eV (actually we stopped our calculation at 11 layers). Other structures are 
thermodynamically unstable. However since the growth is kinetically prevented from forming 
infinite height films, metastable structures with low energies (e.g. finite thickness hR2 films) 



will appear in actual growth. However, the HCP-based structures lie systematically above 
the hR2 ones. By studying the bulk Bi energy with different lateral strains, we confirni that 
the underlying reason that hR2 films are more favorable than HCP films is that Bi favors a 
bilayer structure with strong covalent bonding, (see Appendix Al.) 
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FIG. 5: (Color online) Surface free energy. Dashed-dotted lines stand for monolayer or less. Solid 
and dotted lines stand for HCP and hR2 films on (3x3) cells respectively. Different slopes indicate 
different coverages according to Eq. [3l Stable sequence is indicated by arrows for monolayer or 



more. 



The experimental surface monolayer structure is not clear, as the authors published an 
erratum 10| regarding their 7-layer film being (8x8) rather than the initially claimed (7x7). 
We find both (7x 7) and (8 x 8) surface monolayer are stable from our total energy calculation. 
If the observed surface monolayer is (8x8) rather than (7x7), then a 3 layer hR2 film would 
fit the experimental deposition rate better than the HCP film would (see Appendix A3.). 

Besides the HCP and hR2 (001) films, we also studied the energy of free standing hR2 
(012) films. The 1 and 2 bilayer hR2 (012) films are more stable on Si(lll) than hR2 (001) 
films. However, with the Ni lattice constant rather than Si, the commensurate bilayer (001) 
film has lower energy by 30 meV/atom, and thus the (012) films are not favorable. This is 
also consistent with the experimental observation that no pseudocubic structure appears. 



CONCLUSION 

We study the growth of Bi on Ni(lll) surface using first principles calculations. The 
proposed HCP films pucker under relaxation and are energetically and mechanically unstable 



8 



to adding capping atoms. We find instead that bulk-like (OOl)-oriented hR2 films above the 
surface monolayer are more energetically favorable than HCP films. The hR2 films seem to 
fit with experimental observations (LEED pattern and coverage) equally well as HCP. If our 
model is correct, growth on Ni(lll) might provide a useful synthesis of uniform hR2 bilayers, 



which have been shown to act as two-dimensional topological insulators |l9l . |20|. We hope 
our theoretical work can trigger more interesting work, both theoretical and experimental 
in this subject. 

Besides the surface growth, phenomena of Bi at Ni interfaces also attract attention re- 
cently. Liquid Bi penetrates and segregates at Ni grain boundaries forming bilayer struc- 
tures 2l| in a stable grain boundary phase called a complexion 22|, |23|]. These bilayer 
interfacial structures can possibly explain the long standing puzzle of the liquid metal em- 
brittlement. However, the underlying mechanisms of bilayer segregation and their relation 
with embrittlement have not been revealed at the quantum level. Our study of Bi on Ni 
surfaces serves as a precursor to this interfacial study. In particular, we note that a pair 
of surface monolayer films, one on each surface at a grain boundary, provides an attractive 
model for the observed Bi bilayers. 
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FIG. 6: (Color online) Bi bulk energy of HCP and hR2 structures. Black is the relaxed hR2 
structure. Red is the relaxed HCP structure. Blue is the evenly spaced HCP structure. Points 
(a),(b),(c),(d) correspond to the structures whose energies are plotted in Fig. [71 

APPENDIX 

Al. Lateral strain 

The proposed Bi HCP films have in-plane lattice constant a = 3.7 A (3x3) and 3.8 A ([3- 
112]), while the hR2 films commensurate to Ni(lll) has a = 4.3 A. To explain why the hR2 
films are more favorable than the proposed HCP films we investigate the bonding character 
of Bi by calculating bulk Bi energies with different a. We adopt bulk cells with a (1x1) 
unit cell in the xy plane and 6 atomic layers in the z direction. The a's are fixed while 
lattice constants c's are fully relaxed. For a less than 3.9 A , both the relaxed HCP and hR2 
structures are evenly spaced. For larger a values, Bi in both structures pairs up to bilayer 
structures. Fig. [6] shows Bi HCP and hR2 bulk energies with fixed a. Clearly, the paired 
hR2 bilayer structures with a ^ 4.6 A are more favorable than the proposed HCP structures 
with a ^ 3.5 A. By paring up, the Bi chemical bonds change from metallic to covalent in 
nature. This strongly affects the periods and amplitudes of surface energy oscillation of free 
standing Bi films. 
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A2. Quantum size effect 

In thin metallic films, electrons are confined in the verticle direction. At low temperature, 
the energies of the confined electrons varying with film thickness governs the relative stability 
of the films. This leads to "electronic growth" |2|, a type of QSE. Based on the usual 
quantization rule, the energy oscillation period in a free electron model is half of the Fermi 



wavelength 2J]. In a solid, taking account of the band structure, the actual energy oscillation 



is the superposition of different oscillations at high symmetry points in the surface Brillouin 



zone 251]. For HCP metal, the (001) electron confinement energy can be written as. 



E{N) = Ap sin {2kfNd + 0f ) + Am sin {2kMNd + (f)^) + A^ sin {2kiiNd + 0^) (4) 

where F, M and K are three high symmetry points in the surface Brillouin zone, /c's are the 
Fermi wave vectors, A's characterize the importance of those three points, 0's are the phase 
shift of three oscillations, A^ is the number of layers and d is layer spacing. The resulting 
energy is the superposition of those three oscillations. We compare the band structure 
prediction and the total energy oscillation of Bi films. 

Fig. [7] shows the surface energy oscillation with different a values. Here the surface energy 
is defined as, 

Isurf = ^l^tot - EBiNBi] (5) 

where Efot is the energy of the film, Nsi is the number of Bi atoms in the film, A is the 
surface area and Esi is the linear part of the total energies of the film structures as in Eq. [TJ 
Shown in Fig. El^a), for HCP films with a = 3.5 A, the oscillation is complex in both period 
and amplitude refiecting the superposition of multiple periodicities. Fourier analysis of '-/surf 
yields a period of around 3 layers. From our band structure calculation, the bulk Fermi 
wave vectors at F, M and K are, respectively, 0.597r/d, O.SSn/d and O.lSn/d. Averaging 
over these three frequencies, the resulting oscillation period is 2.9 layers, which agrees well 
with the total energy calculation 'jsurf- This QSE due to electron confinement does govern 
the energy oscillation of the free-standing HCP film with a = 3.5 A. However, the energy 
oscillation amplitude is much smaller than the energy oscillation of Bi HCP films on Ni(lll), 
implying that quantum size effect is not the dominant factor in determining the stability of 
Bi HCP film on Ni(lll). Also unlike Bi on Ni(lll), SOC does aher the relative stability of 
free standing HCP films with a = 3.5 A. 

12 



Fig. [Tl^b) shows the surface energies of HCP films with a equal to 3.9 A where the bulk 
Bi atoms start to pair up. The energy oscillation period is close to 2 layers which is different 
from the band structure prediction of 3.2 layers. Fig. Wic) and El^d) show the energies of 
HCP films with a equal to 4.4 A and hR2 films with a equal to 4.6 A respectively. In both 
cases, the energies show bilayer oscillation with much larger amplitude than the oscillation 
due to this QSE. We conclude that the oscillations of 'jsurf are not due to this QSE for 
a > 3.9 A instead it is due to covalent bonding into bilayers which is a different type of 
QSE. 
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FIG. 7: (Color online) Surface energies of free standing Bi films with different in-plane lattice 
constants a. Black and red curves are using collinear calculation. Blue using noncollinear. 
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A3. Deposition rate 

We analyze Fig. 1 in the support material of [l| to estimate the coverage of the 3x3 film. 
As the deposition rate is constant, the total number of Bi atoms doubles from deposition 
time of 6000s to 12000s. The surface monolayer has fractional area of around 90% at 6000s 
and around 20% at 12000s. Since the fractional area of [3-112] is small, we just ignore the 
small coverage difference between it and the (3x3) film in our estimation. We then get the 
equation, 

0.9 X e„,ono + 0.1 X ^3x3 = -(0.2 X O^ono + 0.8 X ^3^3) (6) 

where 9mono and 6*3x3 are the coverage of surface monolayer and 3x3 film respectively. We 
then get, 

6*3x3 = 2.7 X Omono (7) 

The resulting coverage of the (3x3) film is ^3x3 = 1.4 Bi/Ni if the surface monolayer is (7x7) 
with 9mono = 25/49 = 0.56 and 6*3x3 = 1.1 if the surface monolayer is (8x8) with 9mono = 
25/64 = 0.39. In comparison, the 3 layer (3x3) HCP film has coverage ^3x3 = 12/9 = 1.3 
and the 3 layer (3x3) hR2 film has coverage of ^3x3 = 10/9 = 1.1. The hR2 film thus fits 
with the deposition rate better than the HCP film if the surface monolayer is (8x8), while 
the HCP film fits better if the surface monolayer is (7x7). 
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